Simulating quantum Brownian motion with single trapped ions 
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We study the open system dynamics of a harmonic oscillator coupled with an artificially engineered 
reservoir. We single out the reservoir and system variables governing the passage between Lindblad 
type and non-Lindblad type dynamics of the reduced system's oscillator. We demonstrate the 
existence of conditions under which virtual exchanges of energy between system and reservoir take 
place. We propose to use a single trapped ion coupled to engineered reservoirs in order to simulate 
quantum Brownian motion. 
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I. INTRODUCTION 

The dynamics of closed systems may be calculated ex- 
actly by solving directly the Schrodinger equation. In 
realistic physical conditions, however, the system one is 
interested in is coupled to its surrounding. In this case 
the dynamics of the total closed system can be extremely 
complicated. For this reason, since the very early days 
of quantum mechanics, a huge deal of attention has been 
devoted to the study of the dynamics of open quantum 
systems 0. 

Nowadays the interest in such a broad field has no- 
tably increased mainly for two reasons. On the one hand 
experimental advances in the coherent control of single 
or few atom systems have paved the way to the realiza- 
tion of the first basic elements of quantum computers, 
c-not 2] and phase quanturngates Q. Moreover, the 
first quantum cryptographic [4j and quantum teleporta- 
tion |5| schemes have been experimentally implemented. 
These technological applications rely on the persistence 
of quantum coherence. Thus, understanding decoherence 
and dissipation arising from the unavoidable interaction 
between the system and its surrounding is necessary in 
order to implement real-size quantum computers and 
quantum technologies. On the other hand, one of the 
most debated aspects of quantum mechanics, namely the 
quantum measurement problem, can be interpreted in 
terms of environment induced decoherence 0. Accord- 
ing to this interpretation the emergence of the classical 
world from the quantum world can be seen as a decoher- 
ence process due to the interaction between system and 
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environment |8|. 

A paradigmatic model of the theory of open systems 
is a harmonic oscillator linearly coupled with a reservoir 
modelled as an infinite set of non interacting oscillators. 
Indeed this model is central in many physical contexts, 
e.g., quantum field theo ry [i| , quantum optics [TL [Tot fTl| ■ 
and solid state physics [l2l |. 

In order to describe quantitatively and qualitatively 
how the reservoir affects the system dynamics one needs 
to make some assumptions on its nature and properties. 
Some of these properties, as for example the tempera- 
ture of the reservoir, can be experimentally measured. 
Other parameters, as the reservoir spectral density or 
the system-reservoir coupling, are assumed on the basis 
of physical reasonableness and deduced by the compari- 
son with experimental data. In this sense the model is a 
phcnomcnological one. 

The importance of the damped harmonic oscillator is 
also due to the fact that it is one of the few exactly solv- 
able non-trivial systems. In fact the Heisenberg equa- 
tions of motion for the total system (oscillator plus reser- 
voir) can easily be solved. The solution of the Heisenberg 
equations of motion is indeed the most straightforward 
method for the description of the dynamics of the expec- 
tation values of observables of interest, e.g. the mean 
energy of the system. An example showing the easiness 
and conceptual transparency of the method is given in 



Another way to describe the time evolution of the sys- 
tem is to look at the reduced density matrix which is 
obtained by tracing the density matrix of the total sys- 
tem over the environmental degrees of freedom. This 
procedure is motivated by the fact that, in general, one 
is interested only in the dynamics of the system and not 
in the state of the reservoir. An exact master equation 
for the reduced density matrix can be formulated and 
exactly solved [3 13 13 13 13 13 El 
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During the last decade huge advances in laser cooling 
and trapping experimental techniques have made it pos- 
sible to confine harmonically a single ion and cool it down 
to very low temperatures where purely quantum mani- 
festations begin to play an important role |22|. A single 
laser cooled ion is theoretically equivalent to a particle 
moving in a harmonic potential, whose center of mass 
motion is quantized as a harmonic oscillator. Such a 
system is a unique experimental system since it approx- 
imates very well a closed system |22j. Indeed unwanted 
dissipation, which in this case manifests itself as a heating 
process depopulating the vibrational ground state of the 
ion, is negligible for times much longer than the usual 
times in which experiments take place [51 13- More- 
over, arbitrary states of the ion motion can be prepared 
and coherently manipulated using proper laser pulses 
[25l I2II . Even extremely fragile states as the Schrodinger 
cat states have been realized and detected [27j • Quite re- 
cently, by using multiple ions in a linear Paul trap, exper- 
iments on quantum non-locality have been performed 
and many-particle entangled states have been realized 
|29j | . Furthermore cold trapped ions are the favorite can- 
didates for a physical implementation of quantum com- 
puters 0,0. 

The aim of this paper is to study the interaction of a 
quantum harmonic oscillator with engineered reservoirs. 
In the context of trapped ions it is possible not only to 
engineer experimentally an "artificial" reservoir but also 
to synthesize both its spectral density and the coupling 
with the system oscillator |3Ctl3l| . This makes it possible 
to think of new types of experiments aimed at testing the 
predictions of fundamental models as the one of quantum 
Brownian motion (QBM ) (o r its high T limit: the famous 
Caldeira-Leggett model [l5|)- 

The mean vibrational quantum number of the ion, also 
called heating function [23j , is the central quantity we in- 
vestigate. We study the heating dynamics of the single 
oscillator in correspondence to different reservoir charac- 
teristic parameters. We analyze the influence of the vari- 
ations of the engineered reservoir parameters, e.g., the 
cut-off frequency of the reservoir spectral density, on the 
heating process. The idea of looking at the variation of 
the open system dynamics induced by the changes of the 
relevant reservoir parameters is in fact rather unusual. 
Up to recently, indeed, only dissipation and/or decoher- 
ence due to interaction with the "natural" reservoir were 
studied 0, El El El El HI ■ 

We compare our analytic results for the observable 
quantities of interest, as the heating function, with 
the non-Markovian wave function (NMWF) simulations 
0,13 1 finding a very good agreement. The main result 
of the paper is the experimental proposal for observing 
new features of quantum Brownian motion with single 
trapped ions. We demonstrate that with currently avail- 
able technology a new regime of the open system dynam- 
ics, characterized by virtual phonon exchanges between 
the system and the reservoir, may be explored. 

The paper is organized as follows. In Sec. [H] we in- 



troduce the Master Equation for QBM and its solution 
obtained by means of a superoperatorial approach. In 
Sec. IIIII we study the behavior of the heating function for 
different values of the reservoir parameters. In Sec. llVl 
we review the basic ingredients of the experimental proce- 
dures for engineering artificial amplitude reservoir, as the 
one we study in the paper, and for measuring the heating 
function. In Sec. we describe our experimental pro- 
posal for revealing non-Markovian dynamics of a quan- 
tum Brownian particle, simulated with a single trapped 
ion. Finally in Sec. I VII conclusions arc presented. 



II. EXACT DYNAMICS OF A QUANTUM 
BROWNIAN PARTICLE 



A. Generalized Master Equation 

The dynamics of a harmonic oscillator linearly coupled 
with a quantized reservoir, modelled as an infinite chain 
of quantum harmonic oscillators, may be described ex- 
actly by means of a generalized Master Equation of the 

form [EElEliil 



dpsjt) 
dt 



ih 



= -K%ps(t) A(t)(X^ - n(t)X*P 



>(t)(X 



2\S 



z 7 (t)X s P s p s (t). (1) 



We indicate with X S ( E ) and P s ( s ) the commutator (an- 
ticommutator) position and momentum superoperators 
respectively and with Hq the commutator superoperator 
relative to the system Hamiltonian. The time depen- 
dent coefficients appearing in the Master Equation can 
be written, to the second order in the coupling strength, 
as follows 



A(t) 

7(*) 
U(t) 

r(t) 



k(t) cos(woT)dr, 



h{t) sin(a;oT)dr, 



k(t) sin(o;oT)dr, 



//(r) cos(woT)dr, 



where 



and 



K (r)=a 2 ({E(r),E(0)}), 



p(T)=ia 2 ({E(T),E(0)}), 



(2) 
(3) 
(4) 
(5) 

(6) 
(7) 



are the noise and dissipation kernels, respectively. In 
the previous equations we indicate with a the system- 
reservoir coupling constant, with loq the frequency of the 
system oscillator and with E the generalized reservoir 
position operator. 
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The Master Equation JIJ is local in time, even if 
non-Mar kovian. This feature is typical of all the gen- 
eralized Master Equations derived by using the time- 
convolutionless projection operator technique 0, |34j or 
equivalent approaches such as the superoperatorial one 
presented in l2ij ■ 

The time dependent coefficients appearing in Eq. 
contain all the information about the short time system- 
reservoir correlation. The coefficient r(t) gives rise to a 
time dependent renormalization of the frequency of the 
oscillator. The term proportional to j(t) is a classical 
damping term while the coefficients A(t) and H(t) are 
diffusive terms. 

In what follows we study the time evolution of the 
heating function (n(t)} with n quantum number operator. 
The dynamics of (n(t)) depends only on the diffusion 
coefficient A(i) and on the classical damping coefficient 
7(0 E3- Furthermore, the quantum number operator 
n belongs to a class of observables not influenced by a 
secular approximation which takes the Master Equation 
into the form |H ^ 



dp s (t) A(t)+ 7 (t) 



dt 2 

A (*)~7(*) ro-t 
2 



\2aps{t)a^ - a)ap s (t) - ps{t)a)a\ 
2a*ps{t)a - aa)p s {t) - psi^aa^] 



(8) 



For this reason, in order to calculate the exact time evo- 
lution of the heating function, one can use the solution 
of the approximated Master Equation ©. In the pre- 
vious equation we have introduced the bosonic annihi- 
lation and creation operators a — (X + iP) j \/2 and 
a) = (X — iP) j V2, with X and P dimensionless posi- 
tion and momentum operator. Note that the above Mas- 
ter Equation is of Lindblad type as far as the coefficients 
A(t) ± 7(t) are positive plj . 



B. Analytic solution 



The exact analytic expression for the time evolution of 
the heating function can be obtained from the approxi- 
mated RWA solution. In the RWA the QCF is [3] 



Xt(0 



= e -A r (i)|C|- 



Xo 



,-T(t)/2 p -iu t 



(11) 



with xo QCF of the initial state of the system. The quan- 
tities Ar(t) and T(t) appearing in Eq. Ijllfl are defined 
in terms of the diffusion and dissipation coefficients A(t) 
and 7(4) respectively as follows 



r(t) 



A r (t) 



2 / -/(h) dh, 



-r(t) 



T(«i) 



A(tijdti 



(12) 
(13) 



Eq. I|ll|l shows that the QCF is the product of an ex- 
ponential factor, depending on both the diffusion A(t) 
and the dissipation 7(f) coefficients, and a transformed 
initial QCF. The exponential term accounts for energy 
dissipation and is independent of the initial state of the 
system. Information on the initial state is given by the 
second term of the product, the transformed initial QCF. 

Having in mind Eq. 1)11(1 and using Eq. I)l()[l. one gets 
the following expression for the heating function 

(n(t)) = e - r ^(n(0)> + \ (V r W - l) + A r (t). (14) 

The asymptotic long time behavior of the heating func- 
tion is readily obtained by using the Markovian station- 
ary values for A(f) and j(t). For a thermal reservoir, one 
gets 

(n(t)) = e- rt (n(0)) + n((J ) (l - e"") . (15) 

In the next section we will discuss in detail the dynamics 
of the heating process and we will show the changes in 
the short time dynamics due to the variations of typical 
reservoir parameters such as its temperature and cut-off 
frequency. 



The superoperatorial Master Equation can be ex- 
actly solved by using specific algebraic properties of the 
superoperators [Tflj . The solution for the density matrix 
of the system is derived in terms of the Quantum Charac- 
teristic Function (QCF) Xt(0 at time t, defined through 
the equation 



Ps(t) 



1 

2^ 



Xt(Oe(^-t'«)d% 



(9) 



It is worth noting that one of the advantages of the super- 
operatorial approach is the relative easiness in calculating 
the analytic expression for the mean values of observables 
of interest by using the relation 



(a tm a"> 



d 



(-^) e' 5l " /2 x(?) 



(10) 



5=0 



III. NON-MARKOVIAN DYNAMICS OF 
LINDBLAD AND NON-LINDBLAD TYPE 

In a previous paper we have presented a theory of 
heating for a single trapped ion interacting with a nat- 
ural reservoir able to describe both its short time non- 
Markovian behavior and the asymptotic thermalization 
process |38j . In this paper we focus instead on the case of 
interaction with engineered reservoirs. In the trapped ion 
context, it is experimentally possible to engineer "artifi- 
cial" reservoirs and couple them to the system in a con- 
trolled way. Since the coupling with the natural reservoir 
is negligible for long time intervals (2^, Hj| , this allows to 
test fundamental models of open system dynamics as the 
one for QBM we are interested in. In a sense this extends 
the idea of using trapped ions for simulating the closed 
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dynamics of quantum optical systems |2J, 1391 |40j to the 
possibility of simulating the dynamics of an ubiquitous 
open system as the damped harmonic oscillator. In par- 
ticular, by using the analytic solution, one can look for 
ranges of the relevant parameters of both the reservoir 
and the system in correspondence of which deviations 
from Markovian dissipation become experimentally ob- 
servable. 

In the experiments on artificially engineered ampli- 
tude reservoirs [8ll | the high temperature condition 
Hu> /KT <C 1 is always satisfied. For this reason here 
we concentrate on this regime of the parameters. For 
an Ohmic reservoir spectral density with Lorentz-Drudc 
cut-off 



2w 



TT LO 2 + U> 



2 ' 



(16) 



the dissipation and damping coefficients j(t) and A(t), 
appearing in the Master Equation (JHJl , to second order in 
the coupling constant, take the form 



7 (t)= 



9 2 

a L0o r 



1-e" 



: cos(w i)- re~" ct sin(w t) , (17) 



and 



A(t) = 2a 2 kT- 



1 + r 2 
- (1/r) sin(wot)]} 



{1- 



,— U) c t 



cos(cjoO 



(18) 



with r — uj c /ljo. Equation IjlSJI has been derived in the 
high T limit, while j(t) does not depend on temperature. 
Comparing Eq. (|18|l with Eq. (|17fl . one notices immedi- 
ately that in the high temperature regime, A(t) 3> j(t)- 
Having this in mind it is easy to prove that the heat- 
ing function, given by Eq. I|14|l . may be written in the 
following approximated form 



(»(*)> * A r (t), 



(19) 



where we have assumed that the initial state of the 
ion is its vibrational ground state, as it is actually the 
case at the end of the resolved sideband cooling process 
HH HE E3 • For times much bigger than the reservoir 
correlation time tr — l/u) c the asymptotic behavior of 
Eq. (|19fl is given by Eq. 1)150. This equation gives ev- 
idence for a second characteristic time of the dynam- 
ics, namely the thermalization time tt = 1/r, with 
r = a 2 LUQr 2 1 (r 2 + 1). The thermalization time depends 
both on the coupling strength and on the ratio r — uj c /loq 
between the reservoir cut-off frequency and the system 
oscillator frequency. Usually, when one studies QBM, 
the condition r ^> 1, which corresponds to a natural 
flat reservoir, is assumed. In this case the thermaliza- 
tion time is simply inversely proportional to the coupling 
strength. For an "out of resonance" engineered reservoir 
with r <C 1, tt is notably increased and therefore the 
thermalization process is slowed down. 



A further approximation to the heating function of 
Eq. H19|) can be obtained for times t <C ry: 



(n(t)}~ / A(ti)dt x 
Jo 



- (r 2 -l)[l - e^*cos(w t)]-re^ c *sin(w t)j(21) 

This approximation shows a clear connection between the 
sign of the diffusion coefficient A(i) and the time evolu- 
tion of the heating function before thermalization. The 
diffusion coefficient is indeed the time derivative of the 
heating function. We remind that, since A(t) ^> j(t) for 
the case considered here, whenever A(t) > the Mas- 
ter Equation JHJ is of Lindblad type, whilst the case 
A(t) < corresponds to a non-Lindblad type Master 
Equation. From Eq. i|2U|) one sees immediately that while 
for A(t) > the heating function grows monotonically, 
when A(£) assumes negative values it can decrease and 
present oscillations. 

To better understand such a behavior we study in more 
details the dynamics for three exemplary values of the 
ratio r between the reservoir cut-off frequency and the 
system oscillator frequency: r>l, r = l and r< 1. As 
we have already noticed the first case corresponds to the 
assumption commonly done when dealing with natural 
reservoir while the last case corresponds to an engineered 
"out of resonance" reservoir. 

For r» 1 the diffusion coefficient, given by Eq. I|18fl . 
is positive for all t and r since 



A(£) a 1 - e~" c * cos(w *) > 0. 



(22) 



Therefore the Master Equation is always of Lindblad type 
and the heating function grows monotonically from its 
initial null value. Equation l|2U|) shows that, for times 
t < t r and for r > 1, (n(t)) ~ (a 2 Lu c kT)t 2 , that is the 
initial non-Markovian behavior of the heating function is 
quadratic in time. 

For r = 1, a similar behavior is observed since also in 
this case A(i) is positive at all times. 

Finally, in the case r< 1, Eq. (|18|l shows that, if r is 
sufficiently small, A(t) oscillates acquiring also negative 
values. It is worth noting, however, that the long time 
asymptotic value of A(i) is always positive. Whenever 
the diffusion coefficient is negative, the heating function 
decreases, so the overall heating process is characterized 
by oscillations of the heating function. The decrease in 
the population of the ground state of the system oscilla- 
tor, after an initial increase due to the interaction with 
the high T reservoir, is due to the emission and subse- 
quent reabsorption of the same quantum of energy. Such 
an event is possible since the reservoir correlation time 
tr = 1/lj c is now much longer than the period of os- 
cillation Tg = 1/ujo. We underline that, although the 
Master Equation in this case is not of Lindblad type, it 
conserves the positivity of the reduced density matrix. 
This of course does not contradict the Lindblad theorem 
since the semigroup property is clearly violated for the 
reduced system dynamics [l|. 
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IV. EXPERIMENTAL TECHNIQUES 

This section gives a brief review of the experimental 
procedures for engineering artificial reservoirs and for 
measuring the heating function of the trapped ion. Start- 
ing from the careful analysis of the recent experiments 
presented in this section, we will describe, in the Sec. [V] 
our experimental proposal for simulating QBM with sin- 
gle trapped ions. 



previous equation. Therefore, for high T, the reservoir 
discussed in the paper can be realized experimentally 
by filtering the random field, used in the experiments 
for simulating an infinite-bandwidth reservoir, with a 
Lorentzian shaped low pass filter at frequency uj c . The 
change of the ratio r thus would be accomplished simply 
by changing the low pass filter. 



B. Measurements of the heating function 



A. Engineering Reservoirs 

Let us begin discussing the technique used to engineer 
an artificial reservoir coupled to a single trapped ion. 
Reference |3ll ] presents recent experimental results show- 
ing how to couple a properly engineered reservoir with a 
quantum oscillator, namely the quantized center of mass 
(cm.) motion of the ion. These experiments aim at 
measuring the decoherence of a quantum superposition 
of coherent states and Fock states due to the presence of 
the reservoir. Several types of engineered reservoirs are 
demonstrated, e.g., thermal amplitude reservoirs, phase 
reservoirs, and zero temperature reservoirs. 

A high T amplitude reservoir is obtained by applying 
a random electric field E whose spectrum is centered on 
the axial frequency lu z /2tt — 11.3MHz of oscillation of 
the ion |3l| . The trapped ion motion couples to this field 
due to the net charge q of the ion: Hint = —qx ■ E, with 
x = (X, Y, Z) displacement of the cm. of the ion from its 
equilibrium position. Remembering that E oc 52 . ti(bi + 
&]), with bi and b\ annihilation and creation operators of 
the fluctuating field modes, and that X cx (a + a)) one 
realizes that this coupling is equivalent to the bilinear 
one assumed to derive Eq. QJ. 

The random electric field is applied to the endcap elec- 
trodes through a network of properly arranged low pass 
filters limiting the "natural" environmental noise but al- 
lowing deliberately large applied fields to be effective. 
This type of drive simulates an infinite-bandwidth am- 
plitude reservoir [3l| . It is worth stressing that, for the 
times of duration of the experiment, namely At = 20/is, 
the heating due to the natural reservoir is definitively 
negligible 31]. 

The reservoir considered in our paper is a thermal 
reservoir with spectral distribution given by 



I (a) = J(w)[n(w) + l/2] 



7T Loj. + W 



■coth (w/KT), 



(23) 



where Eq. (fCT)|) has been used. For high T, Eq. 
becomes 



2KT dl 

TT LU 2 + LU 2 



(24) 



The infinite-bandwidth amplitude reservoir realized in 
the experiments corresponds to the case u c — > oo in the 



In this subsection we focus on two experimental meth- 
ods for measuring the heating function of a single trapped 
ion. The first method is based on the asymmetry in the 
sideband motional spectrum of the ion and it has been 
used in |23| for measuring the process of thermalization 
of an ion, initially cooled down to its ground vibrational 
state, due to the interaction with the natural reservoir. 
The same method is used in [26] for measuring the cooling 
dynamics of an ion subjected to sideband cooling lasers. 
The second technique allows to measure the populations 
of the vibrational density matrix of the ion, from which 
the heating function can be obtained. This last method 
has been used in [3l| in the case of interaction with an 
artificial amplitude reservoir. 

For both techniques the first step is the preparation 
of the initial vibrational and electronic ground state, 
1 7i = 0, — ) = \n = 0) |— ), obtained by laser cooling 
and optical pumping to the state |— ). The mean vibra- 
tional number is then measured after fixed delay time 
intervals. During the delay artificial noise, which sim- 
ulates the amplitude reservoir, may be applied. In this 
way, the time evolution of (n(t)) is obtained. 

Let us begin with the first technique. At each fixed 
delay time, the ion is in its electronic ground state |— ) 
and in a certain vibrational state. A laser pulse, tuned to 
a vibrational sideband, is then used to transfer the pop- 
ulation to the upper electronic level |+). After this, by 
means of an electron shelving technique, the electronic 
state of the ion is detected in order to check wether a 
transition to |+) has occurred or not. Repeating this pro- 
cedure one gets the electronic excited state occupation 
probability P + . The amplitude of the blue and red vibra- 
tional sidebands is defined as the probability of making 
a transition |— ) — > |+) due to a laser pulse tuned to the 
blue or red sideband, respectively, and therefore is given 
by P + . This quantity depends on the mean occupation 
number (n). For \n = 0), only the blue sideband can be 
exited while the red one is absent. In general the asym- 
metry in the amplitude of the fcth red {I^ ed ) and blue 
{Iblue) vibrational sidebands, allows to extract (n) |23| : 



red 



l 



blue ■ 



(25) 



A limitation of this method is given by off-resonant 
excitation via the carrier transition. If the driving field 
is tuned to the first lower vibrational sideband, in the 
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resolved sideband condition loq ^> fl, with ft Rabi fre- 
quency of the laser pulse, processes involving off-resonant 
transitions go as (fi/wo) 2 - In order to have a sensible 
measurement of the heating function, the population due 
to off-resonant transitions have to be much smaller than 
the scale over which (n(t)) varies. 

We now focus on the second experimental method for 
measuring the heating function. This method actually 
allows to measure the diagonal elements P n of the vi- 
brational density matrix and it has been used to observe 
their decay due to the interaction with an artificial am- 
plitude reservoir, as the one described in previous subsec- 
tion. For this type of reservoir and in the experimental 
conditions of |2j|, the time evolution oi P n (t) = p nn (t) is 
well approximated by the law 

a) 1 y ( "7* y ( i \ 2n ~ 2j 

Pnn{t) ^ j + 2L, y l + n jt J \l + njt) 

(26) 

where the phenomenological parameters n and 7 are the 
mean reservoir quantum number and the heating rate, re- 
spectively. Equation $ZG$ is valid under the assumptions 
of high T reservoir and for times much smaller than the 
thermalization time, yt <C 1. Such conditions turn out 
to be verified in the experiments described in |3lj . In 
this experimental situation, the Markovian behaviour of 
the heating function, before thermalization takes place, 
is simply given by (n(t)) = fvyt. Note that in the limit 
r» 1, corresponding to an infinite bandwidth Markovian 
reservoir, Eq. (|2UJ| becomes (n(t)) ~ 2a 2 KTt/-K. Com- 
paring this expression of the heating function with the 
one of the experiments one can identify n = KT/loq and 
7 = 2a 2 uio/ir. It is worth underlining, however, that ac- 
cording to Eq. H26f) . only the product fvy may be deduced 
from the experimental data. 

In order to describe our experimental proposal for sim- 
ulating QBM, it is important to look in more detail at 
the procedure used in [31| to obtain the time evolution 
of the populations P n . This allows to deduce from the 
experimental data the characteristic parameters of the 
amplitude reservoir used in the experiments. In |3l| . 
after preparing a superposition of Fock states, the am- 
plitude noise is applied for a fixed time of t = 3/zs af- 
ter which the populations P n are measured. In order 
to measure P n the ion is irradiated with a pair of Ra- 
man beams tuned to the first blue sideband for various 
probe times t p , and P-(t p ) is measured from the fluo- 
rescence signal. From the P-(t p ) data the populations 
P n (i) are finally extracted with a single value decom- 
position [33. In order to observe the time evolution of 
P n one can proceed in two equivalent ways. Either one 
changes the interval of time t during which the amplitude 
noise is applied, or one fixes i and varies the variance of 
the applied noise (V 2 ). In fact, as shown in [3l|], the 



variance of the random noise used to simulate an ampli- 
tude reservoir is 2(V 2 ) — fvyt. Practically, increasing the 
fluctuations of the random electric field applied at the 
trap electrodes is equivalent to an increase in the heating 
function (n) = fvyt. Since P n (t), as given by Eq. 1201) . 
depends only on fvyt, one can obtain the time evolution 
of the populations simply by changing (V 2 ). This is the 
method used in |31| for measuring the populations, so, in 
fact, in the experiment P n {(V 2 )) is recorded. In princi- 
ple, by using the relation fvy = (V 2 ) ft, with t = 3^s, one 
could directly obtain the characteristic parameter of the 
reservoir from the value of the noise voltage applied to 
the trap electrodes. However, unknown geometrical fac- 
tors in converting the voltage to variations in the secular 
frequency prevent a direct comparison. If we indicate 
with (V 2 ) ap pi = c(V 2 ) the fluctuations of the voltage 
applied to the electrodes, with eel, fitting the exper- 
imental data according to the theoretical law given by 
Eq. (|26l) , allows to extrapolate the factor c and therefore 
fvy. It is not difficult to show that, in the experiment 
on the decay (heating) of a Fock state due to interac- 
tion with engineered amplitude reservoir (see Fig. [15] of 
[111), c ~ 10 and hence, for t = 3^s and (V 2 ) = 0.25V 2 , 
fvy ~ 0.84 • 10 7 Hz. Note that, in the experiment, (V 2 ) is 
varied from to 0.3U 2 . 

At this point we are ready to describe our experimental 
proposal for observing the non-Markovian dynamics of 
the heating function and, in general, for simulating the 
dynamics of a quantum Brownian particle. 



V. EXPERIMENTAL PROPOSAL FOR 
SIMULATING QBM 

It is well known that non-Markovian features usually 
occur in the dynamics for times t -C tr = l/u> c . In gen- 
eral, since lo c 3> luq and typically uiq ~ 10 7 Hz for trapped 
ions, this means that deviations from the Markovian dy- 
namics appear for times t <C 0.1/^s. This is the reason 
why the initial quadratic behaviour of the heating func- 
tion is not observed in the experiments, wherein the typ- 
ical time scales go from 1 to lOO^s. 

A way to force non-Markovian features to appear is to 
'detune' the trap frequency from the reservoir spectral 
density. This corresponds, for example, to the case in 
which r = lu c /ujq = 0.1. In this case the reservoir cor- 
relation time is bigger than the period of oscillation of 
the ion and this leads to the oscillatory behaviour of the 
heating function predicted by Eq. I|20|l for r<l. Under 
this condition r c = 1/Lts, and therefore the non-Markovian 
features show up in the time evolution and can be mea- 
sured. Detuning the trap frequency from the reservoir, 
however, decreases the effective coupling between the sys- 
tem and the environment and, for this reason, in order 
to obtain values of the heating function big enough to 
be measured we need to increase either the coupling con- 
stant a 2 , which correspond to an increase in the intensity 
of the voltage applied to the electrodes, or the strength 
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FIG. 1: Time evolution of the heating function for 
2a 2 KT/n = 0.84 ■ 10 9 Hz, lo c = 1MHz and r = 0.1. Solid 
line is the analytical and circles the simulation result. 

of the fluctuations (V' 2 ), which correspond to an increase 
in the effective temperature of the reservoir. 

Let us look in more detail to Eq. (|20|l . For r <C I this 
equation becomes 

(n(t)) ~ 2a KT r 2 {uo c t + \l - e~ Uat cos(w t)l 

TTUJc 

-re~^*sin(w i)} ■ (27) 

In the comparison with the experiment done in previous 
section we have seen that the front factor 2a 2 KT/Tr = 
H7 ~ 0.84 • 10 7 Hz. Increasing of two order of magni- 
tude this front factor, and for r = 0.1, Eq. 127|) predicts 
the behaviour for the heating function shown in Fig. ^ 
We believe that for this range of (n) and of times the 
oscillations of the heating function could be experimen- 
tally measurable. For example, with the first technique 
described in subsection II V Bl using £1 = 10 6 Hz and for 
luq = 10 7 Hz, the ground state population transferred to 
the excited level due to off resonant excitation is of the 
order of 10~ 2 , which is one order of magnitude smaller 
than the variation of the heating function we want to 
measure (see Fig. [JJ. Also the other method described 
in previous subsection seems to be enough accurate to re- 
veal the oscillatory behaviour of the heating function in 
the conditions here examined. However, it is worth not- 
ing that, since the heating function does not depend now 
on the dimensionless variable n/yt, but rather on u> c t, the 
time evolution can be obtained only by varying the dura- 
tion of the time of application of the amplitude reservoir. 
This is not equivalent to a change in the applied voltage 
fluctuations, as it was in the Markovian case discussed in 
subsection IIVBI 

Summarizing, in order to observe the virtual exchanges 
of phonons between the system and the reservoir, lead- 
ing to the oscillations of the heating function, one needs 
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FIG. 2: Time evolution of the heating function for 
2a 2 KT/n = 0.84 • 10 9 Hz, lu c = 1MHz, and r = 10. Solid 
line is the analytical and circles the simulation result. 



to increase of two order of magnitude the coefficient 
2a 2 KT/TT — n-y. This can be done either increasing 
the intensity or increasing the fluctuations of the applied 
noise, or combining an increase in the intensity with an 
increase in the fluctuations. Moreover one needs to use 
a low pass filter for the applied noise, as described in 
subsection II V Al having cut off frequency u c = O.lo^o- 

We now examine briefly the conditions for which the 
quadratic behaviour of the heating function, could be 
observed. We remind that this is the case in which 
r 1 and the time evolution of the density matrix is 
of Lindblad-type. In view of the considerations done 
at the beginning of this section, in order to reveal non- 
Mar kovian dynamics in a time scale of 1 — 100/iS we need 
to have tr > 1/xs, that is, for r = 10, ujq < 0.1MHz. 
This means that one actually needs a 'loose' trap. For 
example, for u>q — 100kHz, with the same applied noise 
used in |2l|), i.e. 77,7 = 0.84 ■ 10 7 Hz, the time evolution of 
the heating function is the one shown in Fig [3 

If one wants to perform a fast measurement of (n), 
e.g. using the first method and assuming il = 10 6 Hz, the 
small value of the trap frequency makes it difficult not 
only to implement one of the two methods for measuring 
the heating function, but also to reach the initial ground 
state since the sidebands are not clearly resolved, and 
therefore resolved sideband cooling technique cannot be 
applied. It is worth stressing, however, that contrarily to 
the case of a natural reservoir, which is always 'in action', 
in the case of an engineered reservoir one can switch off 
the applied noise after a certain delay time and, assuming 
that the effect of the natural reservoir is negligible, mea- 
sure the heating function without the severe requirement 
of big values of £1. If one assumes that after an amplitude 
noise pulse the state of the ion does not change, then it 
is not necessary to perform a fast measurement of (n). 
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Therefore one can work with smaller values of f2, such 
that n/w -C 1. 

Concluding, while measuring the quadratic behaviour 
of the heating function (r ^S> 1) could be a more challeng- 
ing task from the experimental point of view, revealing 
the oscillatory non-Markovian behaviour (r <C 1) appears 
to us in the grasp of the experimentalists, in the condi- 
tions we have analyzed in this section. 

VI. CONCLUSIONS 

In this paper we have studied the dynamics of a single 
harmonic oscillator coupled to a quantum reservoir at 
generic temperature T. In our analysis we have used 
both the analytic solution for the reduced density matrix 
and the NMWF method. 

We have paid special attention to the non-Markovian 
heating dynamics typical of short times. In this regime 
the system time evolution is influenced by correlations 
between the system and the reservoir. For certain val- 
ues of the system and reservoir parameters, virtual ex- 
changes of energy between the system and its environ- 
ment become dominant. These virtual processes strongly 
affect the short time dynamics and are responsible for the 



appearance of oscillations in the heating function (non- 
Lindblad type dynamics). 

Extending the ideas of using trapped ions for simulat- 
ing quantum optical systems, we have proposed to sim- 
ulate QBM with single trapped ions coupled to artifi- 
cial reservoirs. We have carefully analyzed the possibil- 
ity of revealing, by using present technologies, the non- 
Markovian dynamics of a single trapped ion interacting 
with an engineered reservoir and we have underlined the 
conditions under which non-Markovian features become 
observable. 
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